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ABSTRACT 

Most millisecond pulsars with low-mass companions are in systems with either helium-core 
white dwarfs or non-degenerate ("black widow" or "redback") stars. A candidate counterpart to 
PSR J1816-I-4510 was identified by Kaplan et al. (2012) whose properties were suggestive of both 
types of companions although identical to neither. We have assembled optical spectroscopy of the 
candidate companion and confirm that it is part of the binary system with a radial velocity ampli- 



tude of 343 ± 7kms , implying a high pulsar mass, Mp 



1.84 ± 0.11 Mq, and a companion 



mass McSin i = 0.193 ± 0.012 Mq, where i is the inclination of the orbit. The companion appears 
similar to proto-white dwarfs/sdB stars, with a gravity logio(5) — 4.9±0.3, and effective temperature 
16, 000±500 K. The strongest lines in the spectrum are from hydrogen, but numerous lines from helium, 
calcium, silicon, and magnesium are present as well, with implied abundances of roughly ten times 
solar (relative to hydrogen). As such, while from the spectrum the companion to PSR J1816+4510 
is superficially most similar to a low-mass white dwarf, it has much lower gravity, is substantially 
larger, and shows substantial metals. Furthermore, it is able to produce ionized gas eclipses, which 
had previously been seen only for low-mass, non-degenerate companions in redback or black widow 
systems. We discuss the companion in relation to other sources, but find we understand neither its 
nature nor its origins. Thus, the system is interesting for understanding unusual stellar products of 
binary evolution, as well as, independent of its nature, for determining neutron-star masses. 
Subject headings: binaries: eclipsing — pulsars: individual (PSR J1816-I-4510) — stars: atmospheres 
— stars: chemically peculiar — subdwarfs — white dwarfs 



1. INTRODUCTION 

The number and variety of millisecond pulsars (MSPs) 
has expanded dramatically in recent years, due to the 
combination of significant radio surveys (e.g., Bailes et al. 
2011; Keith et al. 2012; Boyles et al. 2012; Lynch et al. 
2012) and the advent of the Fermi satellite (see Ray et al. 
2012 for a recent review). These millisecond pulsars of- 
fer unprecedented laboratories to study particle acceler- 
ation (Abdo et al. 2009), binary evolution (Benvenuto, 
De Vito, & Horvath 2012), and the masses of both white 
dwarfs (Corongiu et al. 2012) and neutron stars (Demor- 
est et al. 2010; Romani et al. 2012). 

Quite a number of discoveries have in fact been the re- 
sult of the confluence of these trends, picking out uniden- 
tified 7-ray sources as priority targets for ongoing ra- 
dio surveys. This technique resulted in the discovery 
of the eclipsing MSP PSR J1816-I-4510 (Kaplan et al. 
2012) as part of the ongoing Green Bank North Celes- 
tial Cap (GBNCC) survey (Stovah et al. 2013, in prep). 
This source has a spin-period of 3.2 ms, an orbital period 
of 8.66 hr, and shows both 7-ray pulsations and radio 
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eclipses for « 10% of the orbit. 

While the radio and 7-ray properties resembled a 
black widow or redback system (Roberts 2011), with a 
non-degenerate companion being ablated by the pulsar's 
wind, the blue colors and lack of variability suggested in- 
stead that the companion to PSR J1816-I-4510 was more 
similar to a helium-core white dwarf (van Kerkwijk et al. 
2005), even though the inferred radius seemed too large 
and the radio eclipses were unprecedented among those 
systems. 

Here we present optical spectroscopy of the companion 
to PSR J1816-h4510 (Section 3.1), taken with the Keck I, 
Palomar 5m, and William Herschel telescopes. We con- 
firm through radial velocity variations that the source 
identified by Kaplan et al. (2012) is in fact the companion 
(Section 3.2), and that, at first blush, it resembles a hot, 
low-gravity white dwarf (Section 3.1). However, detailed 
analysis shows that the spectrum has strong lines of he- 
lium and metals (calcium, silicon, magnesium, and pos- 
sibly iron). We discuss the implications of the spectrum 
and radial velocity measurements in Section 4, along with 
possible origins for the unique features of this system and 
ways they can be exploited to learn more about neutron 
stars and their companions. We conclude in Section 5. 

2. OBSERVATIONS & REDUCTION 

We observed the likely counterpart of PSR J1816+4510 
with a number of instruments over 2012, as summarized 
in Table 1. We primarily report data from the blue sides 
of the spectrographs, although the red sides were ana- 
lyzed as well. Flux calibrators were observed with the 
same setups (CD 153 for the WHT data, Feige 34 for 
the 2012 April P200 data, BD-t-33°2642 for the Keck 
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Fig. 1. — Stars used for velocity references. Left: optical images of the PSR J1816+4510 field from Kaplan et al. (2012). The radio 
position is indicated with the tick-marks, and the uncertainties are dominated by uncertainties in the absolute astrometry of the optical 
data. We also label the two reference stars used for our spectra. Right: spectra of the two reference stars with best-fit models. Star 1 
(with the data from the 2012 March Keck observations) is at the bottom, along with a 6250 K model from Munari et al. (2005, offset for 
clarity). Above them is star 2 (data from 2012 June Keck observations) and a 5750 K model. 



TABLE 1 

Summary of Spectroscopic Observations 



Date 




Phase'' 


Exposure 


Telescope /Instruments^ 


Star^l 


Velocity** 


(TDB BJD) 






(sec) 






(km/s) 


2012-02-25 


06:28 


0.116 


900 


WHT,ISIS-I-QUCAM,300B 


1 


-337 ± 16 




06:41 


0.141 


600 






-317 ±20 


2012-03-23 


14:56 


0.908 


600 


Keck I,LRIS-B, 600/4000 A 


1 


-379 ± 12 




15:07 


0.930 


600 




-396 ± 12 


2012-04-20 


10:39 


0.999 


600 


P200,DBSP,600/4000 A 


1 


-457 ± 13 


2012-05-19 


07:29 


0.989 


360 


P200,DBSP,600/4000 A 


2 


-444 ± 15 




09:47 


0.255 


360 




-102 ± 14 




10:47 


0.371 


360 






118 ± 14 


2012-05-20 


08:40 


0.896 


360 


P200,DBSP,600/4000 A 


2 


-374 ± 14 




10:03 


0.057 


360 




-444 ± 13 




11:17 


0.198 


360 






-200 ± 14 


2012-05-22 


09:28 


0.531 


240 


Keck I,LRIS-B,400/3400A 


2 


251 ± 15 


2012-06-13 


07:31 


0.266 


240 


Keck I,LRIS-B,400/3400A 


2 


-57 ± 14 




07:35 


0.274 


240 




-25 ± 13 


2012-08-18 


09:34 


0.381 


300 


Keck I,LRIS-B, 600/4000 A 


1 


146 ± 13 




09:40 


0.394 


300 




165 ± 13 



^ Time of the middle of the exposure, corrected to the Solar System barycenter. 
^ Orbital phase, based on the ephemeris of Stovall et al. (2013, in prep.). 

The telescope, instrument, and grating/grism used. In most cases we only discuss the results from the blue sides of the spectrographs. 
ISIS is the Intermediate dispersion Spectrograph and Imaging System (ISIS) on the 4.2-m William Herschel Telescope, used with the 
QUCAM electron-multiplying detectors. LRIS is the Low Resolution Imaging Spectrograph (Oke et al. 1995) on the 10-m Keck I telescope. 
DBSP is the Double Spectrograph on the 5-m Palomar 200-inch (Hale) telescope. 

Reference star for each observation, as labeled in Figure 1. 
" Radial velocity for the pulsar, relative to star 1. The uncertainties include a ±12kms~^ systematic uncertainty. 



data and the 2012 May P200 data). The sHt width var- 
ied between 0'.'7 and I'.'S depending on the conditions. 
For all observations of PSR J1816-I-4510 except for the 
last two Keck observations we did not have the sht at the 
parahactic angle as we instead included both the pulsar 
and a reference star on the slit to provide a velocity stan- 
dard (for the 2012 June Keck observations the reference 
star was instead observed right after the pulsar). We 
used two different stars, as listed in Table 1 and shown 
in Figure 1. 



For all of the observations our reduction followed stan- 
dard procedures, although we used MIDAS for the WHT 
data and IRAF for the other data. We subtracted bi- 
ases and overscans then flattened the data using dome 
flats, where the shape of the lamp was removed with a 
high-order polynomial. We then extracted the spectra of 
the pulsar and the reference star using optimal weight- 
ing, subtracting background as measured on both sides of 
the object spectrum. Wavelength calibration was done 
using arc spectra, and fluxes were calibrated using the 
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Fig. 2. — Spectrum of the optical counterpart of 
PSR J1816+4510, based on the 2012 March Keck observa- 
tion. We show the results from both the blue and red sides of 
LRIS. The spectrum is reasonably consistent with the photometry 
(Kaplan et al. 2012 and Kaplan et al. 2013, in prep.), given the 
quality of the flux calibration. Along the top are ticks for the 
different absorption lines seen in the spectrum. Topmost (green) 
are Balmer lines, as expected from a DA white dwarf. Next (red) 
are lines from He I. Finally (black) are low-ionization metal lines 
(mostly Ca I, Ca II, Mg II, Si II). These are shown in more detail 
in Figure 4. 



TABLE 2 

Elemental Abundances 



Species logio(Af/AfH) logio(A^/A^H)soiar [X/U] 



He 0.0 ±0.5 -1.0 1.0 ±0.5 

Na^'.... < -3.0 -5.7 < 2.7 

Mg -3.8 ±0.2 -4.4 0.6 ±0.2 

Si -3.8 ±0.2 -4.5 0.7 ±0.2 

Ca -4.5 ±0.3 -5.6 1.1 ±0.3 

Feb.... -2.8 ±0.3 -4.5 1.7 ±0.3 



Note. — Solar abundances are from Anders &i Grevesse (1989). 



^ Most of the Na seen in the spectrum was interstellar, with sig- 
nificant wavelength shifts compared to the other species. 
^ Detection of Fe consisted of numerous lines near the noise thresh- 
old. Individual lines were not significant, but the ensemble was. 

spectrophotometric standards, correcting for differences 
in airmass using extinction curves appropriate for each 
site. 

3. ANALYSIS 
3.1. The Spectrum of PSR J1816+4510 

To start, we focus on the highest-quaUty spectra which 
are the 2012 March Keck data. We show a set of blue and 
red-side spectra in Figure 2. One immediately notices the 
broad Balmer lines, consistent with our tentative iden- 
tification of the optical counterpart as a low-mass white 
dwarf (Kaplan et al. 2012). However, it is also clear that 
Balmer lines are not the only ones present. We also see 
strong but narrower lines from He I, with at least 20 
lines present between 3000 A and 7000 A, as well as nar- 
row lines from Ca II, Si II, and Mg II, and possibly Fe I 
(see Figure 4 for a more detailed view). 

We fit the data to white-dwarf model atmospheres to 
constrain the atmospheric parameters log]^o(ff) Tcff 
as well as the He and metal abundances. While we obtain 
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Fig. 3. — Results of model fitting to the companion to 
PSR J1816+4510. For each value of the helium number density 
(indicated by the colors, and labeled by the log of the number den- 
sity relative to hydrogen), we show the best-fit values of effective 
temperature and logarithm of the surface gravity, determined for 
aU LRIS 600 Imm^r observations (filled circles), LRIS 4001mm-i 
observations (open squares), and DBSP observations (stars). De- 
pending on the exact model there is considerable scatter and/or 
covariance among the best-fit values. Overall, the model with 
logio(He) = 0.0 had the lowest x^, followed by log]^Q(He) = +0.5, 
then logj^g(He) = —0.5. 

reasonable results, we note that the validity of our mod- 
els has not really been tested for gravities logio(g) < 5. 50 
where our best fits are. Instead those gravities are more 
typical of main-sequence/sdB stars. With this caveat, 
and acknowledging the inhomogeneous data-set, our fit- 
ting procedure was as follows. 

Our model grid covered gravities logiglg) =4-6 (in 
steps of 0.25 dex) and temperatures T^s = 14,000- 
20,000 K (in steps of 250 K). We also included helium, 
with relative abundances logj^pHe/H = —1 to +1 by 
number. For each spectrum we interpolated the model 
spectrum to the observed wavelength grid for a range 
of trial velocities, convolving the model with the pixel 
width and a Gaussian determined by the seeing, and 
truncating by the slit; we explored also convolving by 
the expected radial velocity change during each obser- 
vation (see below) but that did not change the results 
noticeably. We normalized the continuum using some 
line-free regions and a 3"^ order polynomial fit to each 
spectrum. The x^(logio(5)7 '^eff, w) was found by sum- 
ming over the 3750-4455 A region, excluding metal lines, 
and done separately for each value of He/H. 

After optimizing for velocity (see below), the fit was 
covariant among remaining parameters, with the he- 
lium abundance altering not just the line depths but 
also the continuum shape. We find reasonable fits for 
logipHe/H — —0.5 to +0.5: going to log^Q He/H ~ 
— 1 (i.e., Solar) increases by 8670 compared to 
logjQHe/H = 0, for 9372 degrees-of-freedom. For each 
value of the helium abundance, the best-fit locus in 
{\ogiQ{g),Tes) space was also covariant, with a posi- 
tive correlation between the fitted parameters (Figure 3). 
Overall more helium resulted in a lower effective temper- 
ature and surface gravity. Our best-fit was log^o He/H = 
0.0±0.5, logio(5) = 4.9±0.3, and Tes = 16,000±500K; 
we believe these uncertainties to be conservative. We see 
no evidence for variations in effective temperature with 
orbital phase, down to a limit of about ±500 K (limited 
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Fig. 4. — Normalized spectrum of the optical counterpart of PSR J1816+4510, based on the 2012 March Keck observation. We show the 
results from both the blue and red sides of LRIS. The broad Balmer lines are indicated with green ticljs at the bottom. He I lines have 
red ticks at the top, while metal lines have black ticks and are labeled by species. We overplot (green) our best-fit model spectrum. The 
region between 5200 A and 5500 A where the model disagrees slightly with our data is increasingly affected by uncertain calibration due to 
the dichroic. 



by the non-uniformity of the data and the uneven phase 
saniphng) . 

With the background model determined we examined 
the metal lines. We fit for the abundances of the individ- 
ual species listed in Table 2 (helium was fit along with 
the spectrum, as discussed above). Other elements were 
present in our models but they did not produce signif- 
icant hues in the 3700-7000 A range. While we could 
see lines from sodium in the spectrum, they had differ- 
ent velocities from the rest and therefore are presumably 
mostly interstellar, hence we quote only an upper limit. 
For iron, no individual line is significantly detected, but 
the large number of lines makes the presence of iron very 
likely (Figure 4). There may be other species present at 
low levels, and we could also determine a number of up- 
per limits, but given the model uncertainties and our dif- 
ficulty in interpreting just the current data we leave the 
analysis here until the models and data have improved. 
Overall, we find that the metals have abundances con- 
sistent with roughly 10 x solar abundance relative to hy- 
drogen, or roughly solar relative to helium. 

3.2. The Radial Velocity of PSR J1816+4510 

We determined radial velocities of the companion to 
determine its orbit. Unfortunately, in the blue part of 
the spectrum where the majority of the spectral lines are 



there are not many sky lines that can be used as abso- 
lute velocity references. This deficiency is made more se- 
vere by the two-armed nature of the spectrographs used, 
where the red- and blue-side wavelength calibrations are 
separate. Instead, we oriented the slit to include both 
the pulsar's companion and another reference star. As 
discussed above, the observations used two separate ref- 
erence stars (Table 1). 

Using the models of Munari et al. (2005), we fit for the 
spectral type of our reference stars over the wavelength 
range 3700A-5400A, allowing for an overall scaling to 
account for errors in flux calibration, slit losses, etc. Ref- 
erence star #1 (SI) appears to be of spectral type mid- 
to late-F (effective temperature 6250 K), while reference 
star #2 (S2) is a bit cooler, roughly G2 (5750 K), as seen 
in Figure 1. This is reasonably consistent with the ob- 
served colors (from 2MASS [Skrutskie et al. 2006] and 
other archival data), where S2 is slightly redder than SI. 

To determine the radial velocities of the pulsar's com- 
panion and the reference star, we used appropriate model 
atmospheres, minimizing the between the model and 
the data for a range of trial velocities. At each veloc- 
ity we convolved the model as described above, and al- 
lowed for a low-order polynomial difference between the 
model and the data. We restricted the wavelength range 
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Fig. 5. — Radial velocities of tlie optical counterpart of PSR J1816+4510 plotted against orbital phase, repeated twice for clarity. Data 
from WHT/ISIS are squares, P200/DBSP are circles, and Keels I/LRIS are diamonds. We show data using reference star 1 (open symbols) 
and star 2 (filled symbols). The best-fit circular orbit is also shown, with residuals in the bottom panel. The thick portion near phase of 
0.25 is the approximate phase range where the radio pulsar is eclipsed. 



to 3700A-5400A (the maximum wavelength was 4450 A 
for the WHT data) so as to focus on a weh-cahbrated 
part of the data above the highest-order Balmer hnes 
and Balmer break but below the dichroic cutoffs. 

For the pulsar, the fits were good, with reduced typ- 
ically near 1 for the WHT and P200 data. The higher- 
quality Keck data had reduced of « 4, so we scaled 
the velocity uncertainties such that the reduced was 
1. The fits were also good for the reference stars, al- 
though not quite as good as for the pulsar. Instead 
was typically about 1.5 for the P200 and WHT data and 
up to 6 for the Keck data. This may be from a slight 
mismatch in resolution (although we explored changing 
our resolution), but more likely is from having an imper- 
fect model grid that is coarse and that may not reflect 
the true metallicity or rotation. In any case, the best-fit 
velocities did not depend much on the exact model or 
the choices of the other parameters. 

The use of two reference stars complicates our analysis 
slightly, as each star can of course have its own peculiar 
velocity. We therefore obtained two spectra where SI and 
S2 were on the slit together without the pulsar. These 
exposures were right before the observations of the pulsar 
on 2012 August 18 and used the same configuration and 
reduction. We found a consistent radial velocity offset of 
—54 ± 1 kms^i^ for S2 relative to SI. In what follows, we 
apply this offset to determine the velocities of the pulsar 
relative to SI. 

We show the best-fit radial velocities in Figure 5. The 
systemic velocity is relative to SI (in no case were the 
finite durations of the exposures significant compared to 
the orbital period, with corrections of < Ikms^^^). We 
have included a ±12kms~^ systematic velocity uncer- 
tainty in the data to account for limitations in our re- 
duction and differences between the instruments. This 
will also account for mis- alignment of the slit intended 



to cover two objects (an offset of approximately one half 
of the slit width corresponds to a velocity shift of ap- 
proximately lOkms^i^). With that, we were able to fit 
a circular orbit using the phase established by the ra- 
dio ephemeris (Stovall et al.), which has since been con- 
firmed in 7-rays (Kaplan et al. 2012) and is more accurate 
than an optical-only ephemeris. With all of the data we 
find a radial velocity amplitude — 343 ± 7kms~^ for 
the companion, with x^ = 14-5 for 14 degrees of free- 
dom and a systemic velocity of —99 ± 8kms~i^ relative 
to SI. If we only use the P200 data including S2 — a 
somewhat more homogeneous data-set — we find a con- 
sistent Kc = 339 ± llkms"!^ (x^ = 3.5 for 4 degrees of 
freedom). We can also reference all of the velocities to 
the 5577 A sky line. Here we might worry that the line 
is at the extreme edge of our spectra where the wave- 
length solutions are not as reliable, but the result is con- 
sistent: Kc = 347 ±11 kms""'^ and a systemic velocity of 
-108 ± 12 kms^i (x^ = 20.9 for 12 degrees of freedom). 

4. DISCUSSION 

Combining the best-fit radial velocity amplitude of 
the companion, Kc — 343 ± 7kms~^, with the mea- 
sured projected semi-major axis of the pulsar x^sK = 
0.595402 ± 0.000002 s, we can infer the mass-ratio: 



A/, 



PSR 



Mc 



KcPb 
27ra;psRC 



9.54 ±0.21. 



The minimum masses for the pulsar and companion are 
then: MpsrSin^i = 1.84 ± 0.11 Mq and Mc sin^ i = 
0.193 ± 0.012 Mq, where i is the inclination of the orbit. 
Note that, unlike black widow systems that are strongly 
affected by heating such that the photocenter and cen- 
ter of mass are not the same (van Kerkwijk, Breton, & 
Kulkarni 2011; Romani et al. 2012), the minimal mod- 
ulation seen for PSR J1816+4510 (Kaplan et al. 2012, 
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Kaplan ot al. 2013, in prep), suggests that we are mea- 
suring the true ccnter-of-mass motion. 

Combining this mass determination with our 
spectroscopicahy-determined gravity, we can deter- 
mine the radius of the companion and the distance to 
the system. In what follows, we give the values for 
the minimum mass, equivalent to that for an edge-on 
orbit, which is not unreasonable given the eclipses 
(although since they are radio eclipses from ionized 
material, constraints on the inclination are weak). The 
minimum radius is 0.26 ± 0.08 i?©, implying an average 
density of 15gcm~^. To get the distance, we use the 
bolometric corrections appropriate for this effective 
temperature from Tremblay et al. (2011, using the 
0.2 Mq models)^, with Mboi - R = -1.52 ± 0.09 mag; 
while this is for a larger gravity (log^oCS') = 6.6) than 
what we measure, the dependence of the bolometric 
correction on gravity is quite small. We assume an 
extinction Ay = 0.22 mag, which comes from fitting 
our photometry to the same synthetic photometry as 
we used for the bolometric corrections*. Overall we 
find d = 4.5 ± 1.7 kpc using the i?-band photometry 
from Kaplan et al. (2012), where the uncertainty is 
dominated by the uncertainty on the measured gravity 
and temperature {d oc T°g.'^^10~'°Sio(^)/^), and we get 
similar results integrating our model atmosphere di- 
rectly. The distance is almost a factor of two larger than 
the dispersion-measure distance discussed in Kaplan 
et al. (2012), suggesting that, as is found for pulsars 
out of the Galactic plane generally (Roberts 2011), the 
dispersion-measure distance is an underestimate. The 
nominal inferred radius is correspondingly larger, indi- 
cating that this source is even further from equilibrium 
than previously thought with a thermal timescale of 
« IMyr (the radn of 0.2 WDs are < 0.03 i?©), 
although it is still well within the Roche-lobe radius of 
0.56 i?0. The distance also implies a larger 7-ray lumi- 
nosity, Lj = 5 X lO^^ergs"^, roughly equal to the total 
spin-down luminosity, E = 5 x 10"^^ergs~^ (ignoring 
possible corrections for secular acceleration and assum- 
ing the standard moment of inertia / = lO^^gcm^), so 
Lry/E w l.O^Qy. Note that the < E requirement is 

somewhat rough, with observed values of L^/E (based 
on poorly-known dispersion measure distances, and 
assuming beaming into 47rster) ranging from 10% to 
> 200% (Ransom et al. 2011); given the uncertainty in 

our photometric distance we see no problems. 

For the nominal E, we expect « 1.1 x 10^^ erg s~^ cm ^ 
of incident flux at the surface of the companion. This 

compares with crT^ff = 3.7 x lO'^^ ergs~^ cm coming 
from the companion itself, so we would expect variations 
of at most ±100 K coming from the day and night sides 
of the companion (assuming an albedo of 100%). Our 

Sec http:/ /www. astro. umontrcal.ca/-bcrgcroii/CoolingModcls/. 
* We note that Ay = 0.22 mag is somewhat less than we had 
assumed in Kaplan et al. (2012), but it is more consistent with the 
total Galactic extinction inferred from molecular gas. The large 
values of Ay in Kaplan et al. (2012) were primarily the result of 
fitting the optical and GALEX ultraviolet photometry simultane- 
ously. With better data it appears that the models significantly 
overpredict the ultraviolet photometry, and when we fit only the 
optical data we find values of Ay near 0.2 mag (Kaplan et al. 2013, 
in prep). 



spectra cannot yet constrain such a variation, but this 
should be possible in the future with high signal-to-noise 
spectroscopy or photometry. 

Our very low gravity brings up a question of nomen- 
clature. Initially we called the companion a white 
dwarf, and indeed we fit the spectra with white dwarf 
model atmospheres. But this gravity is really more 
consistent with a sdB star or a proto-WD, such as 
GALEX J171708.5-h675712 (Vennes et al. 2011) or 
HD 188112 (Heber et al. 2003). As shown in Figure 6, 
PSR J1816-I-4510 is even more extreme in the gravity axis 
than those sources. It almost certainly is not a standard 
sdB star, as these typically have effective temperatures 
> 25,000 K and masses of ss 0.5 Mq (able to sustain core 
He fusion) compared to 16,000 K and 0.2 here^. It 
may still be contracting and cooling from a recent hy- 
drogen shell flash, as discussed in Kaplan et al. (2012). 
We note that the approximate mass of 0.2 Mq is at the 
limit where He-core WDs are expected to undergo shell 
flashes instead of steady burning (Panel et al. 2007): if 
we could determine a more precise mass and establish 
whether or not unstable flashes occurred, it will help our 
understanding of the evolution of low-mass WDs. 

Another question is the origin of the metals^". While 
planetesimal accretion has been suggested for other low- 
gravity WDs (Figure 6; Gansicke et al. 2012), this seems 
very unlikely for PSR J1816-I-4510, where the relativistic 
wind of the pulsar and the ionized wind of the compan- 
ion will both suppress accretion. The solar Ca/He ratio 
also supports a primordial origin for the metals (cf. Du- 
four et al. 2010). At the other mass extreme, comparison 
with sdB stars also shows differences. For instance, Edel- 
mann et al. (2003) show that helium abundance increases 
with effective temperature, becoming close to solar only 
for Tefi > 30,000 K. Geier (2012) finds no correlation of 
metal abundances with Tcff or the helium abundance, al- 
though they are generally sub-solar. The ionized wind 
that likely causes the radio eclipses may be part of the 
answer here, as it might alter diffuse equilibrium (much 
as inferred for sdB stars winds; Unglaub & Bues 1998) 
and hence the onset of shell burning. It also may remove 
enough angular momentum to create ultra-compact bi- 
naries such as PSR J1719— 1438 and its planetary-mass 
companion (van Haaften et al. 2012). Alternatively, if 
a shell flash occurred recently, this might not only ex- 
plain the low surface gravity but also the abundances, 
as it might have mixed up the interior and the present 
low gravity would delay sedimentation compared to an 
ordinary white dwarf. 

Another interesting comparison is the black widow sys- 
tem PSR J1311-3430 (Romani et al. 2012). Its compan- 
ion has a lower gravity, a much lower mass, and is likely 
less degenerate than the companion of PSR J1816+4510, 
yet has a comparable maximum temperature and a spec- 
trum dominated by lines of He I, with small amounts of 
metals (Ca II, Mg II) visible and no detectable hydro- 

^ While the inclination is unknown so in principle a companion 
mass of f» 0.5 Mq is possible, our constraint on the mass ratio 
would then require a neutron star mass of > 5 Mq 

If interpreted as a white dwarf, the companion would be a 
DABZ star (hydrogen atmosphere with significant helium and met- 
als) , which is an almost unique combination, especially at this grav- 
ity. The companion is in fact more metal-rich than the previous 
record holder, SDSS J073842.56+183509.6 (Dufour et al. 2010). 
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Fig. 6. — Effective temperature vs. log(surface gravity) for low-mass white dwarfs. The WD-WD binaries (largely from the ELM surveys 
of Kilic et al. 2012) are the triangles/squares. Binaries with A stars (such as from van Kerkwijk et al. 2010) are 6-pointed stars, while 
WD-PSR binaries are 5-pointed stars. All sources with high-quality spectra and log^Q(g) < 5.6 have Ca II in their atmospheres (circles); 
those with T^ff > 15,000 K have other species as well (diamonds). The evolutionary tracks for M < 0.25 Mq are from Panel et al. (2007), 
and show only the terminal cooling, ignoring shell flashes. PSR J1816+4510 is labeled, as are other sources of interest: the "proto-WD" 
systems GALEX J171708. 5-1-675712 (Venues et al. 2011) and HD 188112 (Heber et al. 2003), pulsars (van Kerkwijk et al. 2005) and select 
WD-WD binaries (KiUc et al. 2012). SDSS J111215.82-I-111745.0 is also labeled: it is notable for having logio(g) > 6 but still showing 
Ca II and for having detectable pulsations that can be used to help determine the depth of the convective zone (Hermes et al. 2012). We 
also show some black widow/redback systems, as compiled by Breton et al. (2013), with the addition of PSR J1311— 3430 from Romani 
et al. (2012). 



gen. Could PSR J1816+4510 evolve into a system like 
PSR J1311— 3430? Current suggestions for the origin of 
PSR J1311— 3430 involve substantial radiative stripping 
(Benvenuto et al. 2012), rather than unstable hydrogen 
burning, so perhaps not. Indeed, there are indications 
for a dense wind from PSR J1311— 3430 that appears 
driven by irradiation; irradiation appears less important 
for PSR J1816-1-4510, but detailed calculations are nec- 
essary to determine this quantitatively. 

The mass of the neutron star is considerably above the 
canonical value of 1.4 Mq. This is consistent with the ac- 
cretion necessary to have spun it up to its current period 
(Tauris, Langer, & Kramer 2012) and with the masses 
of a number of other pulsars with low-mass He white 
dwarf companions (van Kerkwijk et al. 2005; Corongiu 
et al. 2012), almost to the level of black- widow systems 
(van Kerkwijk et al. 2011; Romani et al. 2012). While 
our current measurement is not as extreme as some (De- 
morest et al. 2010), the brightness of the companion and 
its narrow spectral lines mean that we can potentially 
get precision of < 1% on if c and hence on the mass ra- 
tio without uncertainty about the photocenter (cf. van 
Kerkwijk et al. 2011). At this point our uncertainty in 
the masses will be entirely dominated by uncertainties in 
the inclination. A Shapiro delay measurement is likely 
impossible because of the eclipses, so we may have to 
rely instead on combining spectroscopy with modeling 
the photometric light-curve (Reynolds et al. 2007; Ro- 
mani et al. 2012); we have begun detailed time-resolved 
photometry to search for irradiation or ellipsoidal varia- 
tions (expected to be roughly 2%). Given the masses we 



infer, we can expect the system to merge due to gravi- 
tational radiation in about 1 Gyr (Shapiro & Teukolsky 
1983); the actual merger time might be shorter if energy 
is dissipated through tidal interactions, and the bloated 
nature of the companion could make this more significant 
than otherwise inferred (Burkart et al. 2012). 

A related possibility is that we might be able to deter- 
mine the rotation of the companion, which would help 
understand the history of mass transfer in the system 
(e.g., Pablo et al. 2012). If rotating synchronously with 
the orbit, we would expect the equatorial velocity to be 
« 35kms~^. This high velocity combined with the low 
surface gravity and the presence of narrow He lines means 
that we might separate the pressure broadening respon- 
sible for the Balmer lines from rotation. Initial estimates 
suggest this might be possible, but again it requires con- 
siderably better data than we have right now along with 
improved atmospheric modeling. 

5. CONCLUSIONS 

We have described an initial spectroscopic study of the 
companion to the millisecond pulsar PSR J1816-f4510. 
Our findings largely agree with the initial speculation 
from Kaplan et al. (2012): that the companion is hot and 
large, consistent with neither black widow companions 
nor with standard white dwarfs. However, we now find 
that it is even more distinct, in showing strong lines of 
helium and metals in its spectrum, together with the 
abnormally low gravity and ionized-gas eclipses. 

We are left with a number of questions. (1) What is 
the lifetime of such an object? It cannot be too small 
compared to typical MSP ages, since otherwise it would 
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be too unlikely to find it among the ^ 100 known MSP 
systems. (2) Why arc there strong metals and helium? 
Are they related? Is the cause of metals similar to that 
of normal low-mass WDs, or is it unique to this system? 
(3) What are the detailed metal abundances? Are they 
all the same relative to solar? (4) Why is it eclipsing? 
A temperature of 16,000 K is generally thought to be 
too low for a wind. Is that related to the pulsar (is it 
ablating)? Or is it a more general property of the com- 
panion? Finally, (5) does the presence of the pulsar sig- 
nificantly affect (through its relativistic wind) the evolu- 
tion of the companion? These questions have important 
consequences for our understanding of white dwarf and 
neutron star evolution, as well as the study of interacting 
binary systems. 

Partially based on observations made with the William 
Herschel Telescope operated on the island of La Palma 
by the Isaac Newton Group in the Spanish Observato- 



rio del Roque de los Muchachos of the Instituto de As- 
trofsica de Canarias. Some of the data presented herein 
were obtained at the W. M. Keck Observatory, which is 
operated as a scientific partnership among the Califor- 
nia Institute of Technology, the University of California, 
and NASA; the Observatory was made possible by the 
generous financial support of the W. M. Keck Founda- 
tion. The authors wish to recognize and acknowledge 
the very significant cultural role and reverence that the 
summit of Mauna Kea has always had within the indige- 
nous Hawaiian community. We are most fortunate to 
have the opportunity to conduct observations from this 
mountain. We thank David Levitan, Assaf Horesh, and 
Yi Cao for assistance with observations; research at Cal- 
tech was funded by a grant from the NSF. We thank an 
anonymous referee for helpful suggestions. 
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